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A COMPUTATIONAL MODEL FOR TRANSDERMAL DRUG DELIVERY 

CROSS REFERENCE TO RELATED APPLICATIONS 
This application claims the benefit of Provisional Application No. 60/441,420, 
5 filed January 21, 2003. 

BACKGROUND 

Field of the Invention : . The present invention relates to the mathematical 
modeling of diffusive systems. In particular, the present invention is directed to a 
computational model for simulating non-linear diffusion of multicomponent systems 

10 through heterogeneous domains. 

State of the Art : Systemic delivery of drugs by percutaneous permeation 
(transdermal drug delivery) offers several advantages over immediate release oral 
dosing. Principal among these are the avoidance of first-pass metabolism in the 
gastrointestinal tract and liver, the possibility of providing constant plasma 

15 concentrations of drug, and the opportunity for multi-day dosing. It is also possible to 
reduce the rate and severity of adverse events. Transdermal delivery devices 
("patches") are designed to deliver drug from th^ patch through the intervening skin 
layers to the underlying micro vasculature. Transdermal patches have been marketed 
for almost 20 years, v^ith perhaps 10 diff^rent-^small molecules (molecular weight < 500 

20 amu) now delivered transdermally. / ^ ; 

A major obstacle to the wider application of transdermal patches is the low skin 
permeability of many drugs, which prevents effective dosing fi-om patches of 
acceptable size. In some cases, however, chemical means permit alteration of the 
structure and djoiamics of the skin to such an extent that adequate drug permeability 

25 can be achieved. This is known as permeation enhancement. 
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In order to develop and optimize transdermal patches capable of delivering a 
desired amount of a particular drug over time, it is important to understand the 
transdermal drug delivery process, particularly where permeation enhancers are 
required. Mathematical models that describe the diffusion of drugs in the skin, as well 
5 as in the drug-delivery device have been developed. See, for example, A. J. Lee. J. R. 
King, and T. G. Rogers, IMA J. Math, Appl Med. & Biol. 13, 127 (1996); A. J. Lee, J. 
R. King and S. Hibbard, IMA J. Math. Appl. Med. & Biol. 15, 135 (1998). However, 
the complexity of the models already developed can lead to cumbersome expressions. 
The finite element method has also been used to model the diffusion and partitioning, 

10 but this method has mainly been limited to treatment of linear diffusion and constant 
partitioning. {See, for example, P. J. Missel, Annals Biomed. Eng., 28, 1307 (2000). It 
would be an improvement in the art, therefore, to provide a model describing a two- 
component diffusion process through heterogeneous domains, taking into account the 
partitioning of the compounds between different domains and possibilities of 

15 permeation enhancement. Such a model would permit further understanding of the 
transdermal diffusion of a drug with a permeation enhancer and could particularly 
valuable in the design and optimization of transdermal patches designed to utilize a 
permeation enhancer to achieve desired drug delivery performance characteristics. 



20 SUMMARY OF THE INVENTION 

In one aspect of the present invention, the present invention includes a method 
of modeling multicomponent nonlinear diffusion in heterogeneous media. In particular, 
the finite element method has been employed to model nonlinear diffusion of two 
components through a heterogeneous media, with the diffusivities and partition 

25 coefficients taking on different values in each compartment. The diffusivities and 
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partition coefficients were modeled as being concentration dependent, and the finite 
element formulation presented here is applicable to general fiinctional forms for the 
diffusivities and partition coefficients. The capability to model concentration- 
dependent partitioning of the substances at the boundaries between the vehicle and the 
5 skin, and between the different layers of the skin is demonstrated using the Lagrange 
multiplier method. The Lagrange multiplier method, with the interface flux as Lagrange 
multipliers, allows a robust treatment of nonlinear partition coefficients. 

Li another aspect of the present invention, the model describing the nonlinear 
diffusion of two components through heterogeneous media created using the finite 
10 element method is applied to model transdermal drug delivery systems. In one 
embodiment, a simple two-layer domain is chosen to model a transdermal drug delivery 
system, with assumptions that the diffusivity and partition coefficient of the drug in the 
skin depend linearly on the enhancer concentration. The finite element results 
demonstrate the modeling capabilities and potentials of the finite element formulation. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 provides a schematic representation of heterogeneous domains of 
diffusion linked by interface flux. 

FIG. 2 provides a schematic representation of the geometry of a transdermal 
20 delivery system, 

FIG. 3 illustrates part of a finite element mesh. 

FIG. 4 provides a graph illustrating drug flux through the lower boundary as a 
fimction of drug diffusivity in the skin as calculated by a model according to the 
present invention. 
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FIG. 5 provides a graph illustrating drug flux through the lower boundary as a 
function of the partitioning of the drug from the patch into the skin, as calculated by a 
model according to the present invention. 

FIG. 6 provides a graph comparing the fentanyl flux calculated using a model 
5 according to the present invention to experimentally observed fentanyl flux with and 
without the presence of lauryl pyroglutamate (LP). 

TABLE 1 describes difflisivities and partition coefficients calculated in runs 1 - 
7 in Example 1. The difflisivities are in imits of cm'^ hr'^ and the partition coefficients 
are dimensionless. 

10 TABLE 2 provides the parameters values for the calculations in Figure 7. The 

difflisivities are in units of cm'^ hr"\ and the partition coefficients are dimensionless. 



DETAILED DESCRIPTION OF THE INVENTION 
The present invention includes a finite element formulation for multicomponent 

15 nonlinear diffusion in heterogeneous media. In the context of the present invention, the 
formulation is applicable to multicomponent nonlinear diffusion from a transdermal 
patch into skin, with diffusion and partition coefficients taking on different values in 
each subregion. The effect of permeation enhancers on the diffusion of the drug is 
modeled through both the concentration dependent diffusivity and partition coefficient. 

20 In addition, a Lagrange multiplier method for treating the nonlinear partition coefficient 
is described. Histologically, biochemically and functionally, the skin can be divided 
into the uppermost stratum comeum and underlying layers of the viable epidermis. In 
the present invention, however, we have simplified the problem by treating the stratum 
comeum and viable epidermis as a single homogenous domain. 
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Mathematical Model 



In one aspect, the present invention is directed to a finite element method for 
modeling the diffusion and partitioning of molecules between two compartments 
having different diffusion and partition coefficients for each diffusing species. The 
case where the partition coefficients are constant is treated in a straightforward manner 
by scaling the concentrations in each compartment by their respective partition 
coefficient, in effect defining a new continuous concentration variable for which the 
diffusion equation is solved. This method cannot be used when the partitioning is 
nonlinear, i.e., when the partition coefficients are concentration dependent. The 
difficulty of modeling a nonlinear partition coefficient, i.e., continuous flux and 
discontinuous concentration at the interface is overcome by the method of subdividing 
the domain into two subdomains, each corresponding to a separate compartment, and 
using independently defined Lagrange multipliers as fluxes on the interface to join the 
two subdomains. Therefore, the two subdomains are linked by the common interface 
flux (with opposite signs in each subdomain), and the constraint on the concentrations 
at the interface is enforced by a separate equation. {See, O. C. Zienkiewicz and R. L. 
Taylor, The Finite Element Method, 4'^ ed., vol. 1, McGraw-Hill, 1989, pp. 373-376). 

We consider a diffusion domain Q, which is subdivided into two domains, Q} 
and D^, with the interface denoted Y\ (Figure 1). The variables in each domain i are the 



superscript in the notations indicates the relevant domain, and the subscript denote the 
component (the drug as component 1 and the enhancer as component 2). The general 
two-component diffusion equation in domain fi' is 



two-component concentration vector c 




and the interface flux X = 
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dt 



-V-^'Vc')=0 on 

c' = ^' on 
-DV^=h^ on 
-Z>'c'^ =A on r, 

Similarly, in domain Q^, setting the interface flux as -X, we can write 



(1) 



— -V-^^Vc^)=0 on 
dt 

= g' on 

-D^c^,n on 



(2) 



where />' = 



is the diffiisivity tensor in each domain. 



An additional set of equations describing the concentration discontinuity of each 
species at the interface (i.e., partitioning) is needed to complete the set of equations. 



1 2 1 



^1 _ ^1 

kl kl 



p--p- on r, 



(3) 



10 



15 



where kf is the partition coefficient of component i in domain j. Equations (3) can be 



written as a matrix equation 
on 



1 I 1 2 /X 

—rc' — -C^ = Q 



(4) 



where ^ = 
k' 



-L 0 



is the partition coefficient matrix in each domain. 



The initial concentration vectors are given as cXt = 0) = c o, and the initial interface flux 
is assumed to be zero. 

Since the unknowns solved for involve both the concentrations and the interface 
flux, which is a derivative of the concentration, this becomes a mixed formulation. 
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The diffusion coefficients and partition coefficients in each subdomain are taken 
to be a function of the concentrations at the position and time (x, t), i.e.. 



10 



15 



D' = 



k' 



D-Ac\A) D[Ac\A 
1 



(5) 



^11 ^1 '^2^ 



where c\=c[{x,t\ c\^c^{x,t) ^ and thus gives rise to a nonhnear diffusion and 
partitioning problem. 

The set of equations (1), (2) and (4) along with the initial conditions describe 
our problem fully. Notice that in addition to the concentrations, which are normally the 
only solution to a diffusion problem, we are also solving for the interface flux between 
the two subdomains. This increases the size of the problem, but allows us to treat 
nonlinear partitioning between heterogeneous domains in a straightforward manner. 
Using the method of weighted residuals, the weak form of the equations (1), (2) and (4) 
can be found to be 



n' dt n' r, r; 



= fv2^^rfQ+ fVv'V Vc'^/Q- fv^^;i^ff'+ fv^Vdr = 0 
dt i k A 



(6) 



I U' ) 



where v 



* \ are the weighting functions, and Vc = jy^'j 



The unknowns and X are discretized in the standard way 
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10 



I 0 yv^ 



(7) 



where N and Nx are the shape functions and c' and 1 are the nodal concentrations and 
interface flux. The weighting functions are also discretized similarly. Note that the 
shape functions for concentrations and flux are different. 
Substituting (7) into (6), we get 



?' =v'''J jN''Ndn—+ jB''D'BdClc'+ jN'^N^drl+ jN''h'dr\ = 0 
[n" dt of r, r; J 



dt r, r? 



. U * 0 ft^ J 





0 






0 










0 


0 





where B = 



Since equations (8 must hold for arbitrary v\ v\ and w they can be written 
equivalently in matrix form as 



(9) 



where = jN^Ndn , K' = \B^D'Bda , = \N^N;,dr , K';, = \nJ ^NdT , and 
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Collecting the unknowns together into the vector c = 

equations (9) can be written compactly as 

Mv-h Kc+ / = 0 
with 



, the system of 



(10) 







0 


o' 






0 














0 




0 


, K = 


0 






, v = - 




> and / = < 






0 


0 


0 








0 




0 




0 



Since 1 does not have a time derivative, it is excluded fi"om the vector v . 

The system of equations (10) is in semidiscrete form, in that they are discrete in 
space but continuous in time. The time of simulation 0 < ^ < ^£ is now subdivided into 
equal time steps such that f = nAt, n = 0, rirs- c'„ =c'{t") and 1„ =i(^") are the 
concentration and interface flux, respectively, at time step n. We use the generalized 
trapezoidal rule as the time integration method, which is also called the predictor- 
corrector method. 



(11) 



Again, i is excluded from this predictor-corrector formulation, i.e., the velocity 
at time step «+l is given by 
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Substituting (9) into (8) gives a set of nonlinear algebraic equations in the 
5 vector of unknowns c„^.,, 

0 = ric„,,)=-^M(fn.i-K..)+Jic„,,+ fn.t (12) 

where r is called the residual. 

The problem is now reduced to finding so that the residual vanishes, 
/.^., r(c„^i)=0 . The nonlinear algebraic equations (12) is solved using the robust 
10 Newton-Raphson method. Newton's method is an iterative procedure, where the 
iteration number is indicated by the superscript k, for the k-th iteration, c* = c^^i . The 
subscript «+l indicating the time step number will be omitted in the following for 
simplicity. The starting value for the unknown is chosen to be the predictor for the 
time step, i.e., c° = . The linearization of (12) results in 

15 > A c = r ^* > f-LM+ A:+ ^c* + -^Iac = 0 (13) 

^ dc ^ "^{aAt dc dc ) ^ ^ 

The equations (13) are solved for the increment in the unknown vector Ac, 
which is then added to the previous iterate 

c'^' =c* + Ac (14) 
20 The new solution is then checked for convergence. 
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For the case where the external forces are independent of the concentrations and 
interface fluxes, the specific matrix form of (13) is given by 
1 



aAt 



X,nonlin 



AS' 
Ac' 

a1 



aAt 



aAt 



-K] 



-M'c 



0 



where the extra stiffness matrices iST^^,,,^ and J^T^ arise due to the nonlinearities of 

the diffusivity and partition matrices and therefore depend on the specific mathematical 
formula for D' and k'. 



dK' 



£'* = J B 



B-'^NB\dClc'* 
dc' ) 



X^onlin 



^ SICES'* - 
dc' ~ 



(15) 



All the matrices are assembled element by element in the usual way. Extension 
of the formulation to more than two layers and more than two components is 
straightforward. 

It should be noted that the variational principle cannot be used to derive the 
correct weak form of the equation in this case because of the non-symmetric nature of 
the partitioning constraint. 

If we write the total potential of the system in the standard way as 
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n*(c,^)=n(c>p-(^i,c'-^c=jrfr (16) 

where the total potential of the unconstrained system is 
n(c)=ijc^f— -v.(/)Vc)V^^+ jc^hdr 

and 



8^8 8 



The stationarity of (14) leads to the set of equations 



r' n ^ 



0 



In contrast to the set of equations (8), equations (17) give rise to a symmetric 
tangent matrix. This is unreasonable physically because it amounts to scaling the flux at 
the interface by the partition coefficient, leading to creationor destruction of material at 
the interface between the two compartments. 

The generalized trapezoidal rule used here for time integration is for linear 
problems, unconditionally stable for a > 0.5, and is accurate to second order when a = 
0.5, and accurate to first order otherwise. A problem with using this algorithm for the 
parabolic equation is the oscillation of the solution. Considering the homogeneous 
model equation 

we see that the amplification factor 
= exp(- ^Ath ^tAz£iM 
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is negative if > l/(l-a). Therefore, for any a 9^ 1, the numerical solution can exhibit 
spurious oscillations. For this reason, a = 1 is used in our calculations. Because this 
results in solutions that are only first-order accurate, a small time step is used to 
5 compensate. 

Application of Mathematical Model to Transdermal System 
Finding a permeation enhancer or combination of enhancers that achieves the 
targeted flux of the drug over the time that the patch is applied is often a fundamental 
task in the development of a transdermal formulations. The various factors that 

10 influence the flux of the drug are the geometry of the patch (e.g., area and thickness), 
the diffusivities in each domain, the partition coefficients of the drug and enhancers 
firom the patch into the skin, and the dependence of the diffusivities and partition 
coefficients of each component on the concentrations of each component. In the 
present invention, we have formulated a mathematical model that allows the diffusivity 

15 and/or partition coefficient of the drug to be dependent on the concentrations of the 
drug itself as well as that of the enhancer, and vice versa. In this section, we 
demonstrate the effects of these nonlinearities on the flux of the drug in a simple model 
problem. 

The diffusion domain Q of consideration is composed of two cylindrical layers, 
20 the dermal patch (layer 1) and the skin (layer 2). A side view of the domain, along with 
its dimensions, is shown in Figure 2. There are two diffusing species, the drug 
(component 1, with concentration ci), and the permeation enhancer (component 2, with 
concentration C2). The lower boundary (F^) of the skin layer acts as a sink for both 
components, = 0, modeling the uptake of the drug and enhancer by the 
25 microcirculation system. The radius of the skin is set to be large enough so that the 
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flux of the components across the sides of the skin domain is negligible. Also, it is 
assumed that no drug or permeation enhancer flow occurs through the top or the sides 
of the dermal patch, therefore simple zero flux boundary conditions are specified at all 
other boundaries (Fa). The dimensions of the patch is typical of an actual transdermal 
drug delivery device, and the skin thickness was chosen to match the thickness of skin 
used in the experimental studies, the results of which we want to compare to our 
calculations. 

The formulation as we have stated above is very general, capable of modeling 
complicated nonlinear diffusivities and partition coefficients of a general mathematical 
form (5). In the present invention, rather than attempting to characterize all the 
different modeling possibilities, two aspects are isolated and the focus is drawn on the 
diffusivity and partitioning of the drug in the epidermis. In one embodiment, the 
mathematical model for our example problem is given as 



D' = 



0 
J_ 

0 



0 












0 






1 












0 " 




~ Ac\ + B 


0 






0 




0 




1 


0 




Uc\ + V 


1 

J 




0 


1 

/C22 . 



(18) 



J_ 
0 



In other words, the cross-diffusivities are taken to be zero (i.e., a concentration gradient 
of one component does not give rise to a flux of the other component), and the 
diffusivity and the partition coefficient of the enhancer are assumed be constants in 
both the patch and skin layers as are the diffusivity and the partition coefficient of the 
drug in the patch. But the diffusivity and partition coefficient of the drug in the skin 
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layer are linear functions of the enhancer concentration. Thus, the permeation enhancer 
is assumed to act (i.e., enhance the drug flux) through increasing both the diffusivity 
and partitioning of the drug in the skin layer. This linear form of the functions Dii^ and 
k\\ IS based on careful analysis of experimental drug fluxes for different drug and 
enhancer concentrations^. 

With this model, using (13), the zy-th partition of the element nonlinear stiffiiess 
matrices (kI^^J and {kI^,J are found to be 



0 ^(a^,Xi+^/,2^.,2)" 

0 0 



N,dQ 



J 



(19) 

-7 '■ — vT^v.yv. ^ 



0 -"^NM 



0 0 



Part of the finite element mesh used for this problem is shown in Figure 3. Since 
the problem is axisymmetric, the mesh is 2-dimensional with the use of axisymmetric 
elements. The mesh was generated with the MSC Patran software (available from 
MSG. Software Corporation in Santa Ana, Califomia) and consists of 275 4-node 
elements with 4 integration points for numerical quadrature. The finite element code 
that we used was developed in our research group, and has been used for several years 
for various classes of problems. The time stepping and iterative algorithms are as 
described in the previous section. 



EXAMPLE 1 

The first series of calculations, runs 1 through 4, was performed for different 
values of the parameter A, keeping the values of all other parameters fixed. These 
values are tabulated in Table 1. For all runs, the initial drug concentration in the patch 
is 0.09 g/cm^, initial enhancer concentration in the patch is 0.08 g/cm^, and no drug or 
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enhancer is initially present in the skin. The resulting drug flux through the lower 

boundary (r^), i.e., the flux entering the bloodstream, is shown in Fig. 5. 

For runs 5 through 7, we kept the value of A at 0.5x10"^ cm'^ hr"\ and varied the 

value of the parameter C/. The other parameter values were all fixed as in the previous 

5 runs, and are tabulated in Table 1 . The units for the diffusivity parameters are in cm'^ 

hr'\ and the partition coefficient parameters are dimensionless (ratios of 

concentrations). The initial concentrations are the same as before, and the resulting 

fluxes through the boundary Fg are shown in Figure 6 as a fimction of time. 

The values of the parameters used in the calculations were chosen somewhat 

10 arbitrarily, due to the difficulty of obtaining experimental measurements. The base 
diffusivities of the drug in the patch and epidermis were selected to match the 
experimentally observed flux of fentanyl through a patch/skin system of the same 
dimensions (see Figure 7). The diffusivity of the permeation enhancer was chosen to 
be significantly smaller than that for the drug based on some experimental indications 

15 (the exact value could not be found). The base partition coefficient of 0.15 of the drug 
in the skin layer is an experimentally measured value. The partition coefficient value 
0.1 of the enhancer lauryl pyroglutamate in the skin is based on experimental 
measurements that yielded values ranging fi*om 0.04 to 0.13. The partition coefficient 
of unity of the drug and enhancer in the patch reflects the fact that we are taking the 

20 patch to be the reference domain firom which the drug and enhancer partition into the 
skin. All other parameters were selected somewhat arbitrarily, to demonstrate the 
sensitivity of the calculated flux to the changes in diffusivities and partition coefficients. 
The magnitude of the parameters A (- 10"^) and £/ (~ 1) were chosen so that the terms 
Ac2 and Ucz would be of the same order of magnitude as the base parameters B and F, 

25 respectively. 
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The results of the runs shown in Figures 5 and 6 are as expected, that is, as the 

diffusivity or the partitioning of the drug into the skin increases, the total drug delivered 

by a certain time (the cumulative flux or area under the flux curves) will increase. The 

concentration profile and the flux through the boundary Tg of the permeation enhancer 

5 are the same for runs 1 through 7, since the enhancer diffusion is linear with the same 

parameter values. Notice that for all the runs 1 through 7 the maximum flux occurs at 

the same point in time (around 1 1 hours) due to this, and the shapes of the flux curves 

are very similar. For comparison, the result for the case where the value of A was kept 

at 0.5x10"^ cm^ hrVconc while the value of B was increased to 1.566 xlO'^ cm"'^ hr"^ is 

10 also included in Figure 5. In other words, the base diffusivity of the drug in the skin 
was increased while the effect of the enhancer concentration on the diffusivity was kept 
the same. It is noted that the flux increases with a steeper slope, reaching the maximum 
flux earlier in time, due to the fact that the diffusivity has increased uniformly, and the 
flux curve is on the whole shifted to the left. Compare this with Figure 6, where the 

15 result for the run in which the value of V was kept at 5, while the value of U was 
increased to 0.2 is shown: the maximum flux is reached at the same time and the flux 
curve do not show a discernible difference in shape. This indicates that the functional 
form of the partition coefficient is less important than that of diffusivity in influencing 
the flux of the drug. 

20 

EXAMPLE 2 

The experimental flux data for the drug fentanyl with and without the 
permeation enhancer lauryl pyroglutamate (LP) is presented in Figure 7. The data with 
fentanyl alone had an initial concentration of 0.08 g/cm^ fentanyl in the patch, whereas 
25 the enhanced data initially had 0.09 g/cm^ of fentanyl and 0.08 g/cm^ of LP in the patch. 

- 17- 



EXPRESS MAIL LABEL NO. EV 325 437 665 US filed January 2 1 , 2004 

ALZ5121USANP (ARC3349) William W. VAN OSDOL, et al. 

The calculated fluxes shown along with the data represent the finite element calculation 
results that most closely match the experimental data firom a series of calculations with 
various parameter value combinations. The initial concentrations for the finite element 
calculations were the same as the initial experimental conditions. The values of the 
5 parameters used for the curves shown are tabulated in Table 2. The units here are again 
cm'^ hr"^ for the diffusivity parameters and dimensionless for the partition coefficient 
parameters. The single component calculation is a linear one, with constant 
diffusivities and partition fimctions. The diffusivity and partition coefficient that 
matched the single component experimental results most closely were then used as the 

10 base diffusivity and partition coefficient for the drug in the two-component calculation, 
resulting in a consist parameter set for the two calculations. Based on this set of 
parameters, a predictive calculation of drug flux with a higher initial concentration 
("loading") of the dmg and enhancer in the patch is also shown in Figure 7. This 
calculation was performed with initial patch concentrations of 0.12 g/cm^ and O.lg/cm"^ 

15 for the drug and enhancer, respectively. It is interesting to note that this moderate 
increase in the initial loading for the drug and enhancer (of around 30 and 25 %, 
respectively) increases the peak drug flux almost two-fold. 

Due to the large number of parameters values that can be varied, along with the 
nonlinear character of the problem, we have not attempted to rigorously fit the 

20 experimental data to the finite element results based on our mathematical model. The 
curves shown here represent the results that most closely match the observed data 
among a series of calculations that we have performed. Even at this level, however, 
comparison of the curves reveals some interesting points. The single set of parameters 
of run 9 in Table 2 captures the rise of the flux to its peak value reasonably well, for 

25 both the single component and two-component diffusion. In contrast, the rapid 
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decrease in flux in the tail part of the observed curves are not reproduced as well in 
either case, although the discrepancy is slightly smaller for the multicomponent 
diffusion. As this rapid drop in flux after reaching the peak is present (in fact, even 
more pronounced) in the single component linear diffusion, it seems that our 2-layer 
5 model of the patch and skin with linear Pick's law is too simple to describe the 
transdermal diffusion process fully. The skin, although it is treated as a homogeneous 
continuum domain in this work, actually has a complex structure consisting of the 
stratum comeum, the epidermis and the dermis layers (didn't mention the dermis in the 
intro), all of which present a different environment to the diffusing drug and enhancer 

10 components. In addition, the patch used in the experiments has an occlusive backing, 
causing water retention, which results in swelling of the patch as well as increased 
hydration of the skin covered by the patch. The effect that this increased water content 
has on the diffusion process is currently not well understood, but it is highly probable 
that it would result in changes in the diffusion and partition coefficients, perhaps with 

15 time, as the hydration level increases. 
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